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ABSTRACT 

In the standard galaxy formation scenario plasma clouds with a high thermal energy 
content must exist at high redshifts since the proto-galactic gas is shock heated to 
the virial temperature, and extensive cooling, leading to efficient star formation, must 
await the collapse of massive halos (as indicated by the massive body of evidence, 
referred to as downsizing). Massive plasma clouds are potentially observable through 
the thermal and kinetic Sunyaev-Zel'dovich effects and their free-free emission. We 
find that the detection of substantial numbers of galaxy-scale thermal SZ signals is 
achievable by blind surveys with next generation radio telescope arrays such as EVLA, 
ALMA and SKA. This population is even detectable with the 10% SKA, and wide 
field of view options at high frequency on any of these arrays would greatly increase 
survey speed. An analysis of confusion effects and of the contamination by radio and 
dust emissions shows that the optimal frequencies are those in the range 10-35 GHz. 
Predictions for the redshift distributions of detected sources are also worked out. 

Key words: galaxies: formation - galaxies: high-redshift - instrumentation: inter- 
ferometers 



1 INTRODUCTION 

A satisfactory theory of galaxy formation requires a good 
understanding of the complex physical processes governing 
the collapse of primordial density perturbations and the 
early galaxy evolution. Measurements of the galaxy lumi- 
nosity and stellar-mass functions up to substantial redshifts 
have highlighted that these functions show conspicuous dif- 
ferences with respect to the halo mass functions predicted 
by the cold dark matter (CDM) theory with the "concor- 
dance" cosmological parameters. At the low-mass end, the 
halo mass functions is much steeper than the galaxy lumi- 
nosity function. As discussed by many authors (Larson 1974: 
Dekel & Silk 1986; Cole 1991; White & Frenk 1991; Lacey & 
Silk 1991; Kauffmann et al. 1993; Cole et al. 1994; Somerville 
& Primack 1999; Granato et al. 2001; Benson et al. 2003), 
the relative paucity of low-luminosity galaxies may be at- 
tributed to the quenching of star formation in low-mass ha- 
los by energy injections from supernovae and stellar winds, 
and by photoionization of the pre-galactic gas. This leads to 
the conclusion that efficient star formation must await the 
collapse of massive halos. On the other hand, the above pro- 
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cesses have little effect on very massive halos, which, in the 
absence of additional relieving mechanisms, would convert 
too large fractions of gas into stars, yielding too many bright 
galaxies, with wrong metallicities (Thomas et al. 2002; see 
Benson et al. 2003 and Cirasuolo et al. 2005 for discussions 
of the effect of quenching mechanisms). An effective cure for 
that is the feedback from active nuclei (AGNs), growing at 
the galaxy centers (Granato et al. 2001, 2004; Bower et al. 
2006; Croton et al. 2006fl 

During their very early evolutionary phases, massive 
proto-galaxies are expected to contain large amounts of hot 
gas, but the gas thermal history is obscure. According to 
the standard scenario (Rees & Ostriker 1977; White & Rees 
1978), the proto-galactic gas is shock heated to the virial 
temperature, but this view has been questioned (Katz et 
al. 2002; Binney 2004; Birnboim & Dekel 2003; Keres et 
al. 2005), on the basis of independent approaches: analytic 



^ Note that the AGN feedback invoked by Granato et al. is rad- 
ically different from that advocated by Bower et al. and Croton 
et al.. The former is a property of all AGNs and is attributed to 
a combination of radiation pressure (especially line acceleration) 
and gas pressure. The latter is associated to the radio active phase 
of quasars ('radio mode' feedback). 
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methods, a high-resolution one-dimensional code, smoothed 
particle hydrodynamics simulations. The general conclusion 
is that only a fraction, increasing with halo mass, of the gas 
heats to the virial temperature. The hot gas is further heated 
by supernova explosions and by the AGN feedback, and may 
eventually be pushed out of the halo. Keres et al. (2005) and 
Dekel & Birnboim (2006) find that there is a critical shock 
heating halo mass of ~ 10""-10^^ Mq, above which most 
of the gas is heated to the virial temperature, while most of 
the gas accreted by less massive halos is cooler. 

The large thermal energy content of the hot proto- 
galactic gas in massive halos makes this crucial evolutionary 
phase potentially observable by the next generation of astro- 
nomical instruments through its free-free emission and the 
thermal and kinetic Sunyaev-Zel'dovich effects (Oh 1999; 
Majumdar, Nath, & Chiba 2001; Platania et al. 2002; Oh, 
Cooray, & Kamionkowski 2003; Rosa-Gonzalez et al. 2004; 
De Zotti et al. 2004). In this paper we investigate the de- 
tectability of this proto-galactic gas exploiting an up to date 
model. For the purposes of the present analysis, the adopted 
model can be taken as representative of the most popular 
semi-analytic models (White & Frenk 1991; Kauffmann et 
al. 1993; Cole et al. 1994; Somerville & Primack 1999; Ben- 
son et al. 2003), that all adopt a similar mass function of 
dark matter halos, a cosmological gas to dark matter ratio 
at virialization, and assume that all the gas is heated to the 
virial temperature. 

Even if some single dish telescopes have the required 
theoretical sensitivity, especially at mm and sub-mm wave- 
lengths (e.g. LMT/GTM, GET at 3 mm, Rosa-Gonzalez et 
al. 2004) these continuum observations will be hampered by 
fluctuations in tropospheric emission. Interferometric array 
observations offer a better trade off between angular reso- 
lution, sensitivity and control of systematics (Birkinshaw & 
Lancaster 2005) together with larger fields of view, and allow 
us to work at lower frequencies where the sources of contam- 
inations from backgrounds and foregrounds are lower and 
may be better estimated. For this reason we focus mainly 
on the capabilities of next generation interferometers: the 
Square Kilometer Array (SKA)Q, the Atacama Large Mil- 
limiter Array (ALMAJj, the Expanded Very Large Array 
(EVLAj3, the new 7 mm capability of the Australia Tele- 
scope Compact Array (ATCA). We note that the situation 
for galaxy scale SZ detection is quite different from that for 
cluster SZ detection. Cluster SZ signals are stronger and 
have much larger angular scales than optimum for the inter- 
ferometer arrays and are best observed with single dishes at 
high quality sites (Carlstrom, Holder & Reese 2002). 

The outline of the paper is the following: in §[2]we briefiy 
describe our reference model; in §[3l we present our pre- 
dictions for the counts of proto-galaxies seen through their 
free-free emission and their thermal and kinetic Sunyaev- 
Zel'dovich effects; in §|4] we analyze the potential of the 
next generation interferometers for detecting such signals, 
describe possible survey strategies, and discuss possible con- 
taminating emissions and confusion effects; in §[5l we sum- 
marize our main conclusions. 



^ http://www.skatelescope.org/ 

^ http://www.alma.info/ 
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We adopt a A-CDM cosmology with h — 0.71, fim = 
0.27, Qx = 0.73, Qb = 0.04, as = 0.8, consistent with 
the results from the Wilkinson Microwave Anisotropy Probe 
(WMAP) (Spergel et al. 2006). 



2 OUTLINE OF THE MODEL 

We adopt the semi-analytic model laid out in Granato et al. 
(2004), with the values of the parameters revised by Lapi et 
al. (2006) to satisfy the constraints set by the AGN lumi- 
nosity functions. For the reader's convenience, we summarize 
here its main features. 

The model is built in the framework of the standard 
hierarchical clustering scenario, taking also into account the 
resuhs by Wechsler et al. (2002), and Zhao et al. (2003a; 
2003b), whose simulations have shown that the growth of 
a halo occurs in two different phases: a first regime of fast 
accretion in which the potential well is built up by the sud- 
den mergers of many clumps with comparable masses; and 
a second regime of slow accretion in which mass is added in 
the outskirts of the halo, without affecting the central re- 
gion where the galactic structure resides. This means that 
the halos harboring a massive galaxy, once created even at 
high redshift, are rarely destroyed. At low redshifts they are 
incorporated within groups and clusters of galaxies. Support 
to this view comes from studies of the mass structure of el- 
liptical galaxies, which are found not to show strong signs of 
evolution since redshift 2^1 (Koopmans et al. 2006). The 
halo formation rate at z>1.5, when most massive early-type 
galaxies formed (Renzini 2006), is well approximated by the 
positive term in the cosmic time derivative of the cosmo- 
logical mass function (e.g., Haehnelt & Rees 1993; Sasaki 
1994). 

We confine our analysis to galaxy halo masses between 
M^;" ~ 2.5 X 10" Mq, close to the mass scale at the bound- 
ary between the blue (low mass, late type) and the red (mas- 
sive, early type) galaxy sequences (Dekel & Birnboim 2006) 
and ~ lO"'^ Mq, the observational upper limit to 

halo masses associated to individual galaxies (Cirasuolo et 
al. 2005). 

The complex physics of baryons is described by a set 
of equations summarized in the Appendix of Lapi et al. 
(2006). Briefiy, the model assumes that during or soon af- 
ter the formation of the host dark matter (DM) halo, the 
baryons falling into the newly created potential well are 
shock-heated to the virial temperature. The hot gas is (mod- 
erately) clumpy and cools quickly in the denser central re- 
gions, triggering a strong burst of star formation. The radi- 
ation drag due to starlight acts on the gas clouds, reducing 
their angular momentum. As a consequence, a fraction of 
the cool gas falls into a reservoir around the central super- 
massive black hole (BH), and eventually accretes onto it 
by viscous dissipation, powering the nuclear activity. The 
energy fed back to the gas by supernova (SN) explosions 
and AGN activity regulates the ongoing star formation and 
the BH growth. Eventually, the SN and the AGN feedbacks 
unbind most of the gas from the DM potential well. The 
evolution turns out to be faster in the more massive galax- 
ies, where both the star formation and the BH activity come 
to an end on a shorter timescale, due to the QSO feedback 
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whose kinetic power is proportional, according to the model, 
to M^Q. 

Mao et al. (2007) found that, for the masses and red- 
shifts of interest here, the evolution of the hot (virial tem- 
perature) gas mass, taking into account both heating and 
coohng processes, is well approximated by a simple expo- 
nential law 



Mhot(t) = Mhot(0)e-'/'■=-^ 



(1) 



where Mhot (0) = /cosm M^ii is the gas mass at virialization, 
Mvir being the halo mass and /cosm « 0-18 the mean cos- 
mological baryon to dark matter mass density ratio. The 
evolution timescale tcond can be approximated as 

teond « 4 X 10= ( ) ^ yr . (2) 



1012 Mq 

The model proved to be remarkably successful in account- 
ing for a broad variety of data, including epoch dependent 
luminosity functions and number counts in different bands 
of spheroidal galaxies and of AGNs, the local black hole 
mass function, metal abundances, fundamental plane rela- 
tions and relationships between the black hole mass and 
properties of the host galaxies (Granato et al. 2004; Cira- 
suolo et al. 2005; Silva et al. 2004, 2005; Lapi et al. 2006). 



2.1 The virial collapse 

The virial temperature of a uniform spherically symmetric 
proto-galactic cloud with virial mass Mvir (dark matter plus 
baryons) and mean molecular weight /i — {2X + 3/4Y)^^, 
X and Y being the baryon mass fractions in the form of 
hydrogen and helium (we adopt X=0.75 and Y=0.25, no 
metals) is 

y ^ 1 umpG Mvir 

2 fcs _Rvir 

where rup is the proton mass, G the gravitational constant, 
and ks the Boltzmann constant. The virial radius Jivir is 
given by 

where pvir is the mean matter density within iivir 

Pvir = Pcf2mA(l -f- Z)^, (5) 

Pc = 3H§/{8ttG) being the critical density. For a flat cos- 
mology {Qm + f^A = 1), the virial overdensity A can be ap- 
proximated by (Bryan & Norman 1998; Bullock et al. 2001) 

_ IStt^ + 82a; - 39a;^ 
with uj = n(z) — 1, and 

{1 + zfQm 



(6) 



(7) 



n^(i + 0)3 + nA 

In the redshift range considered here {z ^ 1.5), we have 

2/3 

{l + z)K, (8) 



5 X 10" 



Mvir 



1012 Mq 



SO that for the massive objects (2.5 x 10^^ Mq < Mvir < 
IO^^'^Mq) we are dealing with, the only relevant cooling 
mechanism is free-free emission. 



We assume that, after virialization, the protogalaxy has 
a NFW density profile (Navarro, Prenk & White 1997): 



cx{l + ca;)2 
where x = r/Rvir, 

Mvir 



Ps 



(9) 



(10) 



with 



log(l + c)- c/(l + c) 



jff = 6.8-10-''ne(y ZfnOCTT 



and c = 3 (Zhao et al. 2003b; Cirasuolo et al. 2005). 
2.2 The free- free emission 

The free-free luminosity of the protogalaxy is computed in- 
tegrating over its volume the emissivity given by (Rybicki 
& Lightrnan 1979): 

y'5ff(rvir,i^)- 
■ exp(— ftpi'/fcBT'vir) ergs""*^ cm~^ Hz~^, (11) 

where the sum in the brackets is on all the chemical species 

in the gas (only H and He in our case) Z being the atomic 
number and and rii the number densities of electrons and 
of ions respectively, C is the clumping factor, for which we 
adopt the value (C = 7) given by Lapi et al. (2006), hp \s 
the Planck constant and ^fi(Tvir, v) is the velocity averaged 
Gaunt factor. For the latter we adopted the analytical ap- 
proximation formulae by Itoh et al. (2000) in their range of 
validity. Outside such range we used the formula given by 
Rybicki & Lightman (1979): 

n3/2- 



gs 



V3 
7r 



17.7 + In 



(12) 



The gas density is assumed to be proportional to the mass 
density (pgas = /cosmp). The electron number density is 



n, = ^(X + y/2), 
rur, 



(13) 



rrip being the proton mass. The adopted value of the clump- 
ing factor C is assumed to be constant with radius, as in 
the model. This rather crude approximation stems from our 
ignorance of the complex structure of the gas distribution. 
Finally, the flux scales with mass, redshift and frequency 

as 

/ 1 _l_ 7\ 7/2 / n -r \ 2/3 

Sff=6.6 X 10-^ffff[Tvir,Kl + 2)] (^-) 



Mvir 



4.8 X 10^® cm 



exp 



1012 Mq 

-1.9 X 10"®(!//20GHz ^ 



,2/3 



Jy,(i4) 



'^^ / M (Mvir/1012 Mq) 

where is the luminosity distance (Hogg 1999): 
dL = ^{l + z)l , , "'" _ . (15) 



2.3 The thermal Sunyaev-Zel'dovich effects 

The inverse Compton scattering of the Comic Microwave 
Background (CMB) photons by hot electrons produces a dis- 
tortion of the CMB spectrum, known as Sunyaev-Zel'dovich 
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(SZ) effect (Sunyaev & Zel'dovich 1972). The distortion con- 
sists in an increase of photon energies which implies a de- 
crease of the CMB brightness temperature at low frequencies 
(i/ < 218 GHz for Tcmb = 2.728 K) and an increase at high 
frequencies: 



ATcMB 



= (a;cothx/2 — 4)y 



(16) 



Tcmb 

where x — hpu / {kBTcMB) and y is the comptonization pa- 
rameter 



y 



dlricTe, 



(17) 



o"T being the Thomson cross-section and the electron 
meiss. In the Rayleigh- Jeans region (a:: << 1) eq. (I16|l sim- 
plifies to 



ATc 



Tcmb 

The SZ signal corresponds to an unresolved flux 



Stsz ~ 2 



(fcsTcMB)^ 
(/IPC)2 



where 



a;* exp(a;) 
(exp(a::) — 1)' 



gix)Y 



(xcotha;/2 - 4) 



(18) 



(19) 



(20) 



and Y is the surface integral of the comptonization param- 
eter. Stsz scales with mass, redshift and frequency as 



0.6 X 10 



Mvir 



1012 M, 



4.8 X lO^^'cm 



0.24 



© 



Jy, 



5/3 



(21) 



and may be positive or negative depending on the sign of 
g{x). Here we will quote only positive fluxes, taking the ab- 
solute value of g{x). 

For a virialized cloud with Mvir — 10^^ Mq at 2: = 2, 
which has a virial temperature Te ~ 1.4 x 10^ K, a mean 
electron density rie — 10"'^ cm~'^ and a virial radius of ~ 
106 kpc, the comptonization parameter is ~ 10~^, yielding 
a negative flux of ~ 60nJy at 20 GHz, on an angular scale 
of ~ 10". 



2.4 The kinetic Sunyaev-Zel'dovich effects 

The kinetic Sunyaev-Zel'dovich effect is due to scattering 
of CMB photons by an ionized cloud moving with peculiar 
velocity v. The associated flux density is (Carlstrom, Holder 
& Reese 2002): 



S,sz = 2^-tf^^hix) 
c (hpc)^ 



where 



h{x) 



x*' exp(a;) 
(exp(a;) — 1)^ ' 



(22) 



(23) 



Te is the optical depth, Vp is the line-of-sight component of 
the velocity, and the surface integral of Te is carried out over 
the solid angle of the moving cloud. The resulting CMB spec- 
trum is still Planckian, but shifted towards higher (lower) 




ogS[Jy] 



Figure 1. Comparison of the differential source counts at 20 
GHz of thermal (solid lines) and kinetic (dashed lines) Sunyaev- 
Zel'dovich effects and free- free (dotted line). For SZ effects we 
obviously use the absolute value of the flux. The counts of the 
kinetic SZ effect include both positive and negative signals, and 
are therefore a factor of 2 larger than those given by eq. I|33|l . 
The decline of the counts of Sunyaev-Zel'dovich effects at the 
faint end is due to the adopted lower redshift (2 ^ 1.5) and halo 
mass (Mvir ^ 2.5 x 1O"M0) limits. 



temperatures for negative (positive) velocities (where nega- 
tive means towards the observer). The function h(x) has a 
maximum at = 218 GHz where the thermal effect vanishes. 
Sksz scales as 

Sksz = 5.3 X 10"** (|«|/393 km/s) [(1 -f z)/if ■ 
■ (Afvir/ 10^^ Mq) (4.8 X lO^^cm/di)^ [h{x)/Q.12\ Jy. (24) 

Following Sheth & Diaferio (2001) we model the distribution 
function of galaxy peculiar velocities, P{v), as a Gaussian 
core with a„ — 680(1 + z)~^^'^ kms~i, extending up to Vt — 
1742(1 -|-2)~i''2 kms^i , followed by exponential wings cut off 
at iimax = 3000(1 -I- z)~i''2 kms~i. Normalizing the integral 
of P{v) to unity, we have: 



P{v)dv = 5.86 X 10"''(1 4- z)^^' {dv/kms~^) ■ 

exp[— 0.5(u/cti,)^] for < vt 

2.065 exp(— 2.3|u|/i'o) for vt < \v\ ^ Vn 



(25) 



where vq = 1000(1 + z)~^^^ kms~^ . The adopted scaling 
with redshift is that appropriate in the linear regime, when 
the effect of the cosmological constant can be neglected, as 
is the case in the z range of interest here. 



3 SOURCE COUNTS 

The mean differential number counts per steradian are given 
by: 

where n{L, z) is the comoving epoch-dependent luminosity 
function per unit dlogL, dV/dz the comoving volume per 
unit solid angle: 
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SZt 100GHz 
SZk 20GHz 
SZk 100GHz 



logS[Jy] 



Figure 2. Comparison of the differential source counts at 20 GHz 
and 100 GHz of thermal (thin and thick solid lines respectively) 
and kinetic (thin and thick dashed linos) Sunyaev-Zel'dovich ef- 
fects. 



dV 

dz 



(27) 



According to our model, for any given z the free-free lumi- 
nosity of a proto-spheroid and its thermal SZ signal depend 
only on its virial mass. The luminosity function can then be 
straightforwardly computed integrating the formation rate 
of virialized objects, d^N{MviT,z)/d\ogMvirdt, over the du- 
ration of the ionized phase and multiplying the result by 
dlog Mvir/rflog L. To avoid unnecessary complications we 
keep the free-free luminosity constant at its initial value over 
a time, tion, equal to the minimum between tcond [eq. ©l, 
the expulsion time of the interstellar gas, Atburst , determin- 
ing the end of the star formation burst, and the expansion 
timescale, and zero afterwards. This simplifying assumption 
implies that the evolution of the hot gas mass [eq. ([T}], den- 
sity, clumping factor, and temperature are neglected. It is 
motivated by our expectation that the effect on the free-free 
luminosity, hence on the counts, of the moderate decrease of 
the hot gas mass over the time tion is counterbalanced by an 
increase of the mean gas density and of the clumping factor, 
as a consequence of the shocks associated to supernova ex- 
plosions and to the AGN feedback. Also, having neglected 
the contribution to the counts from the free-free emission at 
t > tion, partly compensates the possible overestimate due 
to having neglected the decrease of the gas mass. In any 
case, a more sophisticated calculation does not appear to 
be warranted since, as discussed in ij |4.3l and 14.41 the free- 
free signal turns out to be too weak to be detectable, being 
overwhelmed by emissions associated to star formation. 

Mao et al. (2007) found that Atburst (yr) can be ap- 
proximated as 



Atburst ~ 4 X 10' 

' 1 



forM.i,>lOi2M0 
Mvir/IO^^M©) forA/vir < lO^^M©' ^ ' 

The mass function of ionized protospheroids at the redshift 
z is then: 



/ d7Vion(Mvir 
\ dMvir 



t(z) 
t(z) — ti, 



dt 



,d^N{M,,,,z) 
dM^ir dt' 



(29) 



The formation rate of protospheroids is well approximated 
by the positive term of the derivative of the Sheth & Tormen 
(1999) mass function, (diV(Mvir, z)/dMvir)sT, (Lapi et al. 
2006): 



d^iV(Mvi„ z) 
dM^irdt 



a5c{z) 

<72(A/viO 
/ diV(Mvir 



2p 



5.{z) a^p{A4.i,) + aPS^,^{z) 



\ dMvir 





dSc(z) 


^ ST 


dt 



(30) 



where a — 0.707, p = 0.3, Sc{z) is the critical overdensity for 
the spherical collapse, (T(Afvir) is the rms amplitude of initial 
density fluctuations smoothed on a scale containing a mass 
Afvir. In turn, the Sheth & Tormen (1999) mass function 
writes 



/ dN{Mv,,,; 

\ dAfvir 



A4"2. 

ST vir 



p dAav 



dlnAfvir 



(31) 



where p is the average comoving density of the universe, 
u = [5ciz)/as{M^i,)]'^, and 



1/2 p-ai//2 



.1/2 



(32) 



with A = 0.322. 

The calculations leading to the counts of the thermal 
SZ "fluxes" are strictly analogous. In the case of the kinetic 
SZ effect we need also to take into account the redshift de- 
pendent distribution of peculiar velocities, and we have 



dAA(5k 



dlogSksz 



In(zi) 



dV 



ln(un 



dln{z)z^ I d\n{v)vP{v) ■ 



dNio4Myir{z,v)] dlogAfvii. 

d log Afvir d log Sksz ' 



(33) 



where Vmin is the velocity yielding a kinetic SZ "flux" 
SkSz from a galaxy with the maximum considered mass 
(Afvir = lO"-^Af0) at redshift z, dNian{M^iT, z,v)/dM^i, 
is the differential mass function of proto-spheroidal galaxies 
with peculiar velocity v and redshift z, producing a kinetic 
SZ flux Sksz- As before, dA'^ion(Afvir, z, v)/dMvii is computed 
integrating the formation rate of virialized objects over the 
duration of the ionized phase. Equation (|26[) gives the num- 
ber of either positive or negative kinetic SZ signals. The 
comparison of the differential source counts at 20 GHz in 
Fig. [T] shows that the thermal Sunyaev-Zel'dovich effect is 
dominant above 10~* Jy. The decline of the SZ counts at 
faint flux levels is due to the adopted lower limits to halo 
masses and redshifts (Afvir ^ 2.5 x 10^^ A'Iq and z ^ 1.5). 
The very steep slope at the bright end comes from the high 
halo mass cutoff. The free-free counts are very low, indicat- 
ing that this emission is very hard to detect in the radio. 

As illustrated by Fig. |2] the SZ fluxes increase with in- 
creasing frequency in the Rayleigh- Jeans region of the Cos- 
mic Microwave Background. 
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Table 1. Main properties of next generation interferometers. The 
maximum baseline has been calculated considering that the an- 
gular size, for the galaxies in the intervals of mass and redshift 
we are considering, ranges from 5" to 35", and requiring a ratio 
of 5 between amplitude and noise on the visibilities. 10% SKA 
has the same properties as SKA, but the number of baselines is 
1.25 X 10^ 



The angular resolution of an array of antennas is given 





FULL-SKA 


ALMA 


ATCA 


EVLA 


Frequency (GHz) 


10-20 


100 


35-50 


35 


Bandwidth{GHz) 


4 


4x2 


2x2 


8 


Antenna diam. (m) 


12 m 


12 


22 


25 


Efficiency 


0.8 


0.8 


0.8 


0.8 


T.ys (K) 


50 


50 


60-80 


75 


No. of polariz. 


2 


2 


2 


2 


Min. baseline (m) 


15 


15 


30.6 


30 


Max. baseline (km) 


1.4-0.7 


0.2 


0.4-0.3 


0.4 


No. of baselines 


1.25 X 10^ 


700 


10 


350 



4 PERSPECTIVES FOR SEARCHES OF 

IONIZED PROTO-SPHEROIDAL CLOUDS 

4.1 Next generation mm-wave interferometers 

In Table [T] we have collected some of the main properties of 
next generation radio interferometers working at few cm to 
mm wavelengths. 

The Australia Telescope Compact Array (ATCA) is a 6 
22m-dish array. The technical parameters we use here refer 
to the recently completed upgrade to the 7 mm receivers and 
the increase of the bandwidth from the present 2 x 128MHz 
up to 4 GHz (CABBS). The band ranges between 30 to 50 
GHz, with Tsys increasing from 60 K up to 80 K at the top 
end of the band. The system will be fully operational by 
2008. 

The Atacama Large MiUimiter Array (ALMA) is a 50 
12 m antenna array. The lower frequency band with higher 
priority ranges between 84 and 116 GHz, close to the max- 
imum amplitude in flux of the negative signal of thermal 
Sunyaev-Zel'dovich effect. The array will be operational by 
2012. 

The Expanded Very Large Array (EVLA) is an im- 
provement of the sensitivity, frequency coverage, and res- 
olution of the existing VLA. When completed, after 2013, 
it will use the 27 25m dishes of VLA working in the fre- 
quency range 1-50 GHz with 8 GHz bandwidth per polari- 
sation available in the frequency bands 18-26.5, 26.5-40, and 
40-50 GHz. 

The Square Kilometer Array (SKA) is a titanic project 
for an interferometer whose main technical specification is to 
have at least one square kilometer of detecting area in the 
core region. The highest frequency band should span the 
range 16-25 GHz. Several designs are under consideration. 
The parameters we use refer to the small parabolic dishes 
version, which is the only high frequency design being con- 
sidered. The telescope is expected to be fully operational 
after 2020, but a "10% SKA" is expected to be operating 
as early as 2015. Phased array feeds in the focal plane are 
being considered for the lower frequency receivers. If such 
systems were implemented at the higher frequencies they 
would increase the field of view and hence the survey speed 
by factors of up to 50. 



by 



1.02- 



B 



(34) 



where B is the maximum distance between two antennas. 
The field of view normally corresponds to the Half Power 
Beam Width (HPBW) of an antenna 

HPBW = 1.02-^ (35) 

where D is the diameter of the antenna dish. For a Gaussian 
beam the field of view (FOV) is 

Phased array feedhorns add a multiplying factor to this re- 
lation, increasing by the same factor the sky area covered in 
a single pointing. The noise level in an image is given by 



BJ-sys 



At] y t Nb 



(37) 



where Tsya is the system temperature, A is the antenna sur- 
face area, rj is the system (dish and receiver) efficiency, t is 
the integration time, Nhasc is the number of baselines short 
enough to have full sensitivity to observe objects with size 
between 5" and 35", Av is the bandwidth and npoi is the 
number of polarizations. Considering the NFW profile for 
densities enlarges the range of full sensitivity baselines in 
the visibility space improving the resolution without losing 
too much in sensitivity. We made our calculations using a 
reasonably conservative configuration. 

Reference values of the quantities used in the calcu- 
lations for the instruments mentioned above are given in 
Table [U 

Assuming as detection level of an object emitting flux 



Slim the ratio Siim/o-in 



5, for a given telescope the 



integration time required for each pointing can be obtained 
by inverting eq. (|37|l . The number of pointings necessary to 
cover a sky area As is 



As/ FOV. 



(38) 



If the integral counts of sources scale as S''^ , the number of 
sources detected in a given area scales as t'^''^. For a given 
flux, the number of detections is proportional to the sur- 
veyed area, i.e. to t. Thus, to maximize the number of de- 
tections in a given observing time we need to go deeper if 
P > 2 and to survey a larger area if /3 < 2. The number 
of sources detected above a given fiux limit, S'um, within 
a telescope FOV, Vpov, is straightforwardly derived from 
the source counts. The number of such pointings necessary 
to detect A^s sources is rip = Ai's/7Vpov(Siim) and the cor- 
responding surveyed area is As = FOV. The predicted 
integral counts of thermal and kinetic SZ effect for several 
frequencies, covered by the radio interferometers mentioned 
above, are shown in Figs. |3] and IH respectively. The scale on 
the right-hand side of these figures gives the corresponding 
area containing 100 protospheroids. 

The time necessary to reach the wanted Sum with 
S/N = ^ in a. single pointing, tp, is obtained from eq. (|37|) . 
and the total observing time for detecting A^a sources (ex- 



cluding the slew time) is obviously tprip 

for Ns - 



In Figs. [5] and [6] 

100 as a function 
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Figure 3. Integral counts and sky area required to detect the 
thermal SZ effect of 100 protospheroids (right-hand scale) as a 
function of the absolute value of the "flux" at 20, 35 and 100 
GHz. 
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Figure 4. Same as Fig. |3] but for the kinetic SZ effect. As in 
Fig-in the counts include both positive and negative signals; for 
the latter, S is obviously the absolute value of the flux. 



of the absolute value of the thermal and kinetic SZ limit- 
ing flux for the 4 instruments in Table [T] at the frequencies 
specified in the inset. The curves have minima at the val- 
ues of 5iim corresponding to the fastest survey capable of 
detecting the wanted number of sources. Clearly, it will be 
very time-consuming to detect 100 protospheroids with the 
EVLA, and unrealistic with the ATCA. On the other hand, 
since the 7 mm upgrade of ATCA will be operational already 
in 2008, it will be possible to exploit it to get the first test 
of the present predictions, and possibly to achieve the first 
detection of an SZ signal from a proto-spheroidal galaxy. 

The SKA large effective collecting area allows the de- 
tection of thermal SZ signals of 100 protospheroidal galaxies 
at 20(10) GHz in 1(11) minutes with 7(7) pointings reaching 
S'lim = 10-^-^(10-'^-^) Jy in a 0.46(0.47) deg^ area. The 10% 
SKA requires 100 times more time than the full SKA but 
is still faster than EVLA or ALMA. If phased array feeds 
were available at the higher frequencies they would improve 
these surveying times by a factor of up to 50. 



Figure 5. Total survey time for ALMA, SKA, EVLA and ATCA 
to detect 100 protospheroids in thermal SZ at the frequencies 
specified in the inset. 
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Figure 6. Same as figure [5] but for the kinetic SZ effect. 



4.2 Redshift distributions 

The redshift distributions of thermal and kinetic SZ effects 
are illustrated, for 3 values of Siim, in Figs. [7]and[8l They 
are both relatively fiat, as the fast decrease with increasing 
z of the density of massive (i.e. SZ bright) halos is partially 
compensated by the brightening of SZ signals [eqs. (|2ip and 
(|24|) ]. Such brightening is stronger for the thermal than for 
the kinetic SZ. A consequence of such brightening is that 
the range of halo masses yielding signals above a given limit 
shrinks with decreasing redshift, as the minimum detectable 
halo mass increases. The upper limit on masses of galactic 
halos then translates in a lower limit to the redshift distri- 
bution for bright Sum- 



4.3 Contaminant emissions 

The adopted model envisages that the plasma halo has the 
same size as the dark matter halo, i.e. of order of hun- 
dreds kpc. In the central region (with size of order of 10 
kpc), the gas cools rapidly and forms stars. Bressan, Silva, 
& Granato (2002) obtained a relationship between the star 
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Figure 7. Redshift distribution (in bins of width 5z = 0.1) of 
thermal Sunyaev-Zcl'dovich effects at 20 GHz for the flux hmits 
specified in the inset. 




Figure 8. Same as in Fig 
Zel'dovich effect. 



[7] but for the kinetic Sunyaev- 



formation rate (SFR) and the radio luminosity at 8.4 GHz, 
Ls(8.4GHz): 



Ls (8.4GHz) ~ 3.6 x lO' 



SFR 



ergs 



(39) 



Mo/yr 

This relationship is in good agreement with the estimate by 
Carilli (2001) while the equations in Condon (1992) imply 
a radio luminosity about a factor of 2 lower, at fixed SFR. 
The Granato et al. (2004) model gives the SFR as a function 
of galactic age, tgai, for any value of the the halo mass and 
of the virialization redshift (see, e.g.. Fig. 1 of Mao et al. 
2007). We must, however, take into account that eq. (|39} has 
been derived using a Salpeter (1955) Initial Mass Function 
(IMF). For the IMF used by Granato et al. (2004), the radio 
luminosity associated to a given SFR is higher by a factor 
of 1.6 (Bressan, personal communication). The coefficient in 
eq. (|39|l was therefore increased by this factor. 

The mean rest frame 8.4 GHz luminosity at given Mvir 
and 2vir was then obtained using the corresponding SFR av- 
eraged over tgai in the renormalized eq. (|39p . Extrapolations 
in frequency have been obtained using, as a template, the 



Figure 9. Comparison of the fiux associated with star formation 
(solid lines) with the thermal SZ "flux" (dashed lines) at 20GHz 
(left panel) and at 35GHz (right panel), as a function of the viri- 
alization redshift for four values of the virial mass (log(Mvir)= 
11.5, 12., 12.5, 13.2, from bottom to top). 



fit to the Arp 220 continuum spectrum obtained by Bressan 
et al. (2002; solid line in their Fig. 2). Using the continuum 
spectrum of M 82 (solid line in Fig. 1 of Bressan et al.), the 
other standard starburst template, we get essentially iden- 
tical results. In Fig. |9] we compare the flux associated with 
star formation with the thermal SZ "flux" [eq. (|21|l ]. as a 
function of the virialization redshift, for several values of 
Mvir. For a given halo mass, the ratio of the thermal SZ to 
the contaminating signal increases with frequency (and with 
redshift) as far as the contamination is due to radio emission 
associated to star formation (we do not consider here nuclear 
radio emission, which occurs in <10 per cent of galaxies). 
However, already at 20 GHz the thermal dust emission be- 
comes important for the highest redshift sources. Such emis- 
sion is more steeply increasing with frequency than the SZ 
signal, even in the Rayleigh- Jeans region of the CMB, and 
rapidly overwhelms it at >100GHz. The SZ/contamination 
ratio increases with increasing halo mass; therefore the SZ 
detection is easier for the more massive halos. Thus in the 
range 10-35 GHz the thermal SZ is expected to dominate 
over the contaminating signal at least for the most massive 
objects. 

It must be noted that the star forming regions are con- 
centrated in the core of the spheroids, on angular scales of 
the order or less than 1", for the redshifts considered here. 
Long (>3km at 35 GHz for full sensitivity) baselines obser- 
vation with high sensitivity may be able to resolve the star 
forming region positive signal and subtract it from the im- 
age. To achieve this purpose a good sampling of the shortest 
spacings on the uv plane is necessary together with a good 
sampling of the largest ones: with the latter it might be 
possible to reconstruct the contaminated profile, subtract it 
from the former and produce an uncontaminated SZ profile. 
Again, the SKA at high frequencies seems to be the optimal 
instrument. 
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4.4 Confusion effects 

Further constraints to the detection of SZ effects are set 
by confusion fluctuations. Fomalont et al. (2002) have de- 
termined the 8.4 GHz source counts down to 7.5 /xJy. For 
>S'8.4GHz;$l rnjy they are well described by: 

N{> S) ~ 1.65 X lO^^S'"^ " arcsec"^ (40) 

with S in /ijy. The spectral index distribution peaks at q ~ 
0.75 (5oc i^-°). 

For all but one (SKA 10 GHz) of the considered surveys 
the "optimal" depth for detecting 100 sources corresponds 
to 8.4 GHz flux densities within the range covered by Fo- 
malont et al. (2002), so that the confusion fluctuations are 
dominated by sources obeying eq. (|40p . We then have: 

aLf ^ 0.2 (-^) ^^sr' MJy' (41) 

\8 GHz/ 100 arcsec^ 

where Sd, in /^Jy, is the detection limit and ui is the solid an- 
gle subtended by the SZ signal. Equation (|4H) can be rewrit- 
ten as 

<Tconf " V8.4 GUz) VlOO arcsec^y ^ ' 

yielding a 5aconf detection limit of ~ 4/iJy at 10 GHz and 
of ~ 2.3 fiJy at 20 GHz. For the "optimal" survey depths at 
higher frequencies S'd/cTconf ^ 5, implying that they are not 
affected by confusion noise due to radio sources. 

On the other hand, as noted above, at high frequencies 
the redshifted dust emission from distant star-forming galax- 
ies becomes increasingly important (De Zotti et al. 2005). 
To estimate their contribution to the confusion noise, we 
have used once again the model by Granato et al. (2004), 
with the dust emission spectra revised to yield 850^ counts 
consistent with the results by Coppin et al. (2006), and com- 
plemented by the phenomenological estimates by Silva et al. 
(2004, 2005) of the counts of sources other than high-z proto- 
spheroids (see Negrello et al. (2007) for further details). We 
find, for a typical solid angle lo — lOOarcsec^, 5(Tconf flux 
limits due to these sources of 3, 55, and 190 /^Jy at 20, 35, 
and 50 GHz, respectively. Thus at 20 GHz we have signifl- 
cant contributions to the confusion noise both from the radio 
and from the dust emission; the overall 5(Tconf detection limit 
is — 4 /xjy. At 10 GHz the contribution of dusty galax- 
ies to the confusion noise is negligible, while at 100 GHz 
the confusion limit is as high as 2 mjy, implying that the 
detection of the galactic-scale SZ effect is hopeless at mm 
wavelengths. 

Although high-2: luminous star-forming galaxies may 
be highly clustered (Blain et al. 2005; Farrah et al. 2006; 
Magliocchetti et al. 2007), the clustering contribution to 
fluctuations is negligible on the small scales of interest here 
(De Zotti et al. 1996), and can safely be neglected. 



5 CONCLUSIONS 

In the standard scenario for galaxy formation, the proto- 
galactic gas is shock heated to the virial temperature. The 
observational evidences that massive star formation activity 
must await the collapse of large halos, a phenomenon re- 
ferred to as downsizing, suggest that proto-galaxies with a 
high thermal energy content existed at high redshifts. Such 



objects are potentially observable through the thermal and 
kinetic Sunyaev-Zel'dovich effects and their free-free emis- 
sion. The detection of this phase of galaxy evolution would 
shed light on the physical processes that govern the collapse 
of primordial density perturbations on galactic scales and 
on the history of the baryon content of galaxies. 

As for the latter issue, the standard scenario, adopted 
here, envisages that the baryon to dark matter mass ratio 
at virialization has the cosmic value, i.e. is about an order 
of magnitude higher than in present day galaxies. Measure- 
ments of the SZ effect will provide a direct test of this as 
yet unproven assumption, and will constrain the epoch when 
most of the initial baryons are swept out of the galaxies. 

As mentioned in §[T] almost all semi-analytic models for 
galaxy formation adopt halo mass functions directly derived 
or broadly consistent with the results of N-body simulations, 
and it is commonly assumed that the gas is shock heated to 
the virial temperature of the halo. They therefore entail pre- 
dictions on counts of SZ effects similar to those presented 
here. On the other hand, the thermal history of the gas 
is governed by a complex interplay of many astrophysical 
processes, including gas cooling, star formation, feedback 
from supernovae and active nuclei, shocks. As mentioned in 
§[T1 recent investigations have highlighted that a substan- 
tial fraction of the gas in galaxies may not be heated to 
the virial temperature. Also, it is plausible that the AGN 
feedback transiently heats the gas to temperatures substan- 
tially above the virial value, thus yielding SZ signals ex- 
ceeding those considered here. The gas thermal history may 
therefore be substantially different from that envisaged by 
semi-analytic models, and the SZ observations may provide 
unique information on it. 

We have presented a quantitative investigation of the 
counts of SZ and free-free signals in the framework of the 
Granato et al. (2004) model, that successfully accounts 
for the wealth of data on the cosmological evolution of 
spheroidal galaxies and of AGNs (Granato et al. 2004; Cira- 
suolo et al. 2005; Silva et al. 2005; Lapi et al. 2006). 

We flnd that the detection of substantial numbers of 
galaxy-scale thermal SZ signals is achievable by blind sur- 
veys with next generation radio interferometers. Since the 
protogalaxy thermal energy content increases, for given halo 
mass, with the virialization redshift, the SZ "fluxes" increase 
rather strongly with z, especially for the thermal SZ effect, 
partially compensating for the rapid decrease of the density 
of massive halos with increasing redshift. The redshift dis- 
tributions of thermal SZ sources are thus expected to have 
substantial tails up to high z. 

There are however important observational constraints 
that need to be taken into account. The contamination by 
radio and dust emissions associated to the star formation 
activity depends on mass and redshift of the objects, but 
is expected to be stronger than the SZ signal at very low 
and very high frequencies. We conclude that the optimal 
frequency range for detecting the SZ signal is from 10 to 35 
GHz, where such signal dominates over the contamination 
at least for the most massive objects. It must be noted how- 
ever that contaminating emissions have typical scales of the 
order of those of the stellar distributions, i.e. < 1" at the 
redshifts of interest here (see Fomalont et al. 2006), while 
the SZ effects show up on the scale of the dark matter halo, 
which is typically ten times larger. Therefore arcsec resolu- 
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tion images, such as those that will be provided by the SKA, 
will allow to reconstruct the uncontaminated SZ signal. 

The coexistence of the hot plasma halo, responsible for 
the SZ signal, with dust emission implies that the scenario 
presented in this paper may be tested by means of pointed 
observations of high-2 luminous star-forming galaxies de- 
tected by (sub)-mm surveys. 

Confusion noise is a very serious limiting factor at mm 
wavelengths. Contributions to confusion come on one side 
from radio sources and on the other side from dusty galax- 
ies. At 10 GHz only radio sources matter; a modest extrap- 
olation of the 8.4 GHz fijy counts by Fomalont et al. (2002) 
gives a 5(Tconf detection limit Sd — 4 ^3y, for a SZ signal sub- 
tending a typical solid angle of 100 arcsec^. Fluctuations due 
to dust emission from high- 2; luminous star-forming galaxies 
may start becoming important already at 20 GHz; at this 
frequency, quadratically summing them with those due to 
radio sources we find again Sd — 4/iJy, for the same solid 
angle. On the other hand, the high resolution of the SKA 
and of the EVLA will allow us to effectively detect and sub- 
tract out confusing sources, thus substantially decreasing 
the confusion effects. Beating confusion will be particularly 
important for searches of the weaker kinetic SZ signal. 
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